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We study the distribution function of the three dimensional wormlike chain with a fixed orientation 
of one chain end using the exact representation of the distribution function in terms of the Green's 
function of the quantum rigid rotator in a homogeneous external field. The transverse Id distribution 
function of the free chain end displays a bimodal shape in the intermediate range of the chain lengths 
(1.3Lp, 3.5-Lp). We present also analytical results for short and long chains, which are in complete 
agreement with the results of previous studies obtained using different methods. 
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I. INTRODUCTION 



Polymers with contour length L much larger than the 
persistence length L p , which is the correlation length for 
the tangent-tangent correlation function along the poly- 
mer and is a quantitative measure of the polymer stiff- 
ness, are flexible and are described by using the tools of 
quantum mechanics and quantum field theory If 
the chain length decreases, the chain stiffness becomes an 
important factor. Many polymer molecules have internal 
stiffness and cannot be modeled by the model of flexi- 
ble polymers developed by Edwards The standard 
coarse-graining model of a wormlikc polymer was pro- 
posed by Kratky and Porod [7| • The essential ingredients 
of this model are the penalty for the bending energy and 
the local inextensibility. The latter makes the treatment 
of the model much more difficult. There have been a sub- 
stantial number of studies of the Kratky-Porod model in 
the last half century Q-Jl^l (and citations therein). In re- 
cent years there has been increasing interest in the theo- 
retical description of semiflexible polymers [2lj]-[34j]. The 
reason for this interest is due to potential applications 
in biology (35j (and citations therein) and in research on 
semicrystalline polymers (36j |. 

It was found in the recent numerical work by Lattanzi 
et al. 37], and studied analytically in 38] within the 
effective medium approach, that the transverse distribu- 
tion function of a polymer embedded in two-dimensional 
space possesses a bimodal shape for short polymers, 
which is considered to be a manifestation of the semi- 
flexibility. The bimodal shape for the related distribu- 
tion function of the 2d polymer was also found in re- 
cent exact calculations by Spakowitz and Wang (39|. In 
this paper we study the transverse distribution function 
G(to,x = 0,y,N) of the three dimensional wormlike 
chain with a fixed orientation to of one polymer end using 
the exact representation of the distribution function in 
terms of the matrix element of the Green's function of the 
quantum rigid rotator in a homogeneous external field 
|33| . The exact solution of the Green's function made it 
possible to compute the quantities such as the structure 
factor, the end-to-end distribution function, etc. practi- 



cally exact in the definite range of parameters [29| , [33[ • 
Our practically exact calculations of the transverse dis- 
tribution function of the 3d wormlike chain demonstrate 
that it possesses the bimodal shape in the intermediate 
range of the chain lengths (1.3L p , 3.5L p ). In addition, 
we present analytical results for short and long worm- 
like chain based on the exact formula (JTJ) , which are in 
complete agreement with the previous results obtained 
in different ways [l9[ (WKB method for short polymer), 
[H[ (perturbation theory for large chain). 

The paper is organized as follows. Section |TT] intro- 
duces to the formalism and to analytical considerations 
for short and large polymers. Section [TTTl contains results 
of the numerical computation of the distribution function 
for polymers with different number of monomers. 



II. ANALYTICAL TREATMENT 

The Fourier-Laplace transform of the distribution 
function of the free end of the wormlike chain with a 
fixed orientation to = dr(s)/ds\ s _ Q of the second end is 
expressed, according to [331 ] , in a compact form through 
the matrix elements of the Green's function of the quan- 
tum rigid rotator in a homogeneous external field P s (/c,p) 
as 



G(t o ,k,p)=J2(0\P S (k,p)\l)V2rTTPi(t o n), (1) 

1=0 

where n = k//c, and P s (k,p) is defined by 



P s {k,p) = (1 + ikDM s )~ 1 D, 



(2) 



with D and M s being the infinite order square matrices 
given by 



D 



Sli> 



11' 



Mfr = wi8 hV+1 + w l+1 5 l+ltl >, (3) 



and wi = \/l 2 / (4Z 2 — 1). The matrix D is related to 
the energy eigenvalues of the free rigid rotator, while 
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M s gives the matrix elements of the homogeneous ex- 
ternal field. Since P s (k,p) is the infinite order matrix, 
a truncation is necessary in the performing calculations. 
The truncation of the infinite order matrix of the Green's 
function by the n x rt-order matrix contains all moments 
of the end-to-end chain distance, and describes the first 
2n—2 moments exactly. The transverse distribution func- 
tion we consider, G(to,a; = Q,y,N), is obtained from 
G(to, x, y, z, N), which is determined by Eqs. CD"©, in- 
tegrating it over the z-coordinate, and imposing the con- 
dition that the free end of the chain stays in the x — 
plane. As a result we obtain 



G(t o ,x = 0,y,N) 



dzG(t 0l x = 0, y, z, N) 



kdkG(t ,k\ kz=0 ,N)J (ky), 



(4) 



where Jo (ky) is the Bessel function of the first kind [4(| • 
Taking the z-axis to be in the direction of to yields ton = 
0, so that the arguments of the Legendre polynomials in 
Eq. (fT]) become zero, and consequently only even I will 
contribute to the distribution function (j4]). 



A. Short wormlike chain 

We now will consider the expansion of around the 
rod limit N — ► 0, which corresponds to the expansion 
of G(t,k,p) in inverse powers of p. To derive such an 
expansion, we write D in the equivalent form as 



D = D {} 



1 



-ED 



(5) 



with D a Ll , = p- l 8iy and E u , = 1(1 + 1)/2<5 M *. Further 

we introduce the notation (1 + ik,DM s ) 1 = yoy with yo 
and y defined by 



y Q = (1 + ikD°M s 



E 

y = l + ik—DM s y y. (6) 
P 



The iteration of D and y results in the desired expansion 
of y and consequently of G(to,k,p) in inverse powers of 
p, which corresponds to an expansion of G(to,k, ./V) in 
powers of N. The leading order term in the short chain 
expansion is obtained by replacing D by D° in Eq. (fT]) 
as 



where cos 9 — to n is the angle between the tangent to 
and the wave vector k. The connection of Go(to,k,p) 
with the plane wave expansion is due to the fact that the 
Kratky-Porod chain becomes a stiff rod in the limit of 
small N. We have checked the equivalency between the 
plane wave expansion JS]) and the distribution function 
(|7|) term by term expanding ([7]) in series in powers of the 
wave vector k. The arc length TV is equivalent for stiff 
rod in units under consideration to the chain end-to-end 
distance r. In the r space the plane wave ([5]) corresponds 
to the stiff rod distribution function 



G(t , r, N)\ N ^ = S(x)6{y)6{z - N). 



(9) 



The iteration of D in ([5]) and y in ^ generates an 
expansion of G(to, k,p) in inverse powers of p. The cor- 
rections to the plane wave to order 1/p 3 are obtained 
as 



Gi(t 0l k,p)=£ 



,/.,-// : ' • //,,(-) ir 

P VP. 



E 



ikyo-D°M s y D 
P 



VzTTlPi (t n) 



01 



The above procedure yields for I = the short chain 
expansion of the distribution function of the free Kratky- 
Porod chain, which was studied recently in (34[. Unfor- 
tunately, we did not succeed yet in analytical evaluation 
of Gi(to, k, N). Such computation would be an inter- 
esting alternative to the treatment of the short limit of 
the wormlike chain by Yamakawa and Fujii [l9j within 
the WKB method. Nevertheless, following the consider- 
ation in [34j we succeeded in computing the anisotropic 
moments ((rto)") for small N 



<(rt o r> = N n 1 



n n(5n — 

-N + — 

2 24 



f) 



-N J 



(10) 



The first-order correction coincide with that obtained in 
[l9| using the WKB method, while the second-order cor- 
rection is to our knowledge new. The higher-order terms 
in pop can be established in a straightforward way us- 
ing the present method. Note that the computation of 
((rt )") does not require the knowledge of the full distri- 
bution function G(to, r, N). 



B. Large chain 



G (t ,k,p) = J2 [yoD°] ol V2r+TP l (t Q n). (7) 



The latter coincides with the expansion of the plane wave 

El 



-ikr cos 



(8) 



In studying the end-to-end distribution function for 
large iV we utilize the following procedure. We expand 
first the expression (0\P s (k,p)\l) in powers of k 2 . The 
structure of this expansion for I — is presented in Table 
|TJ The subseries in powers of k 2 in the qth column are 
denoted by T l Thus we have 



G*(M = V2rTT(o\P s (k, P )\i) = Y^T'c 
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TABLE I: 



g ■ 


1 


k° : 


1 




P 2 


k 2 : 


- w l 




IT 


fc 4 : 






— — 


k 6 : 






fc 8 : 


™1 



2 

-L — tyf 

p 

c 4 I 22 or- 4 -r 2 2 



3p 2 



- ■ j24 M.'^ — 12(1 W ^ it- 2 — f - ! 2 "'i "'^ ~ 2 " ; 1 

Sip 2 

-117iu? tu|to|+270ujf +9iu 4 

27p 3 



TABLE II: 



G»(fc,p) 



J fc 2 (45p+4fc 2 ) fc 2 (4725 p 3 +90 fc 4 p+980 fc 2 p 2 +2 fc B ) 

fc 2 /3+P + 15(fc 2 +3p) 2 525(fc 2 +3p) a 
3ik _ l/5ifcp(4fc 2 +45p) 

( t 2 + 3p) 2 
k 2 1 fc 2 (fc 4 +45fc 2 p+280p 2 ) 

J (fc 2 +3p) 2 

! ifc 5 (65fc 2 +294p) 
2025 (fc2 +3p ) 2 •- 

1 508 fc s +2250 fc 6 p 
259875 ( fe 2 +3p ) 2 •■• 



3fc 2 +9p ~<~ 

ik 3 
15fe 2 +45p 

_2_ fc 4 

525 fc 2 +3p 



+ 



+ 



The series T l q with small values I and g possess a simple 
structure and can be summed up. For example Tg and 
are given by 



T[' = 3 



E 



(-i) r 



Wl 



2 i n 2 m — 2 



i 



fc 2 /3+p' 



J 2 



2 m 7,2 m 



m l + 



W 2 fc 2 



oo / \m-l 

E B 

m— 1 

fc 2 (45p + 4fc 2 
15 (fc 2 + 3p) 2 



While T° corresponds to the distribution function of the 
Gaussian chain, T° give the qth correction to the Gaus- 
sian distribution. The inspection of the series for 7f and 
T 2 X shows that they are expressed by T® and as 

7f = -iklf, 
rpi _ Zip 



However, it seems that there is no general recursion re- 

z,p) for 
0) and 



lation for T l q . The results of computations of Gi(k,p) for 



I = 1, 2, 3, 4 by taking into account g = 1,2,3 (Z 
q = 1, 2 (i ^ 0) are summarized in Table HT1 

The inverse Laplace- Fourier transform of G/ (fc, p) given 
in the table yields the expansion of the end-to-end dis- 



tribution function G(R,Q,t) to order 0(t 2 ) as 



G(R, 9, t) = (4tt)- 



2 7rt 



5 

3441 i? 6 



i? 2 33 i? 4 
I 2 " ~ 40 ~W 



3/2 

exp 
6799 i? 4 



3i?< 
"~2T 



1600 i 4 
1089 R 8 329 i? 2 79 
1400 t 5 ' 3200 I 5 " ~ 240 ~W ~ 640 i 2 
3i? 153 i? 3 99 i? 5 25 R, 



.on 100 rt, »» n zo n. . . 
+ ( 2T + 1^^~ 80^~16^ )Pl(c ° Se) 



+( 



li? 2 



2 t 2 

3 i? 3 



961 i? 4 
560 IF 



33 i? 6 67 R 2 
80 ~ 60 T 3 " 
9 i? 4 



)P 2 (cos6) 



1400 i 4 



- r P 4 (cose) 



(11) 



where t — N/2, R = r/2, and is the angle between 
R and to. The latter is in accordance with the result 
by Gobush et. al. [H[ derived in a different way. The 
expansion of G(R, 0, t) for large N can be extended in a 
straightforward way to include higher-order corrections. 



III. NUMERICAL RESULTS 

The computation of the distribution function of the 
polymer with the fixed orientation of one end is per- 
formed by truncating the infinite order matrices in ((T|) 
with the finite ones, and by taking into account the fi- 
nite number of terms in the summation over The in- 
verse Laplace transform of dXJ) is carried out with Maple. 
The results of the calculation of the distribution function 
G(t ,x = 0,y,N), using the truncations with 20 x 20 
order matrices and restricting the summation over the 
quantum number I at I = 8, are given in Fig. [T] The 
results show that the distribution function possesses the 
bimodal shape at intermediate chain lengths within the 
interval N £ 1.3, ...,3.5. We also find that the distribu- 
tion function becomes Gaussian for very short and very 
long chains. At the onset there is a tiny maximum at 
y = 0, which we interpret as remnant of the Gaussian 
behavior of short chains. The maximum at y = for 
3d chain is a rather small effect, which is difficult to be 
explained in a qualitative way. 

We now will discuss qualitatively the origin of the bi- 
modal behavior of the projected distribution function of 
the free end of the wormlike chain. The very short worm- 
like chain behaves similar to a weakly bending stiff rod, so 
that the distribution function of the free end is Gaussian 
with the maximum at y — 0. The typical conformation 
of the chain in this regime looks like a bending rod with 
constant sign of curvature along the chain. For larger 
contour lengths the curvature fluctuations are small and 
are still controlled by the bending energy, however with 
varying sign of curvatures along the chain. The typical 
conformation of the chain can be imagined as undulations 
along the average conformation of the polymer. The pro- 
jected distribution function of the free end in this regime 
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FIG. 1: Normalized distribution function G(to,x — Q,y,N) 
for various chain lengths, computed with 20 x 20 matrices. 
The insets show the distribution function at the onset of bi- 
modality, and in the region of its disappearance. 



is expected to be roughly uniform within some range of 
y. We expect that the inhomogeneities of curvature fluc- 
tuations in the vicinity of the clamped end are the reason 
for the maximum at y ^ 0. The larger curvatures in the 
vicinity of the fixed end result in larger displacements 
y of the free end, and therefore contribute preferentially 
to the maximum at y 7^ 0. With further increase of N 
the conformations correspond to undulations around the 
average conformation of the chain, which is now a me- 
andering line. Fluctuations become now less controlled 
by bending energy, which results in weakening of the bi- 
modality. Since the difference between 2d and 3d chain 
is assumed to be marginal for short chains on the pro- 
jected distribution function, we will compare the onset 
of the bimodality in both cases. Because the transver- 
sal displacement is measured in both cases in units of 
the contour length, we have to recompute the number of 
segments for 3d and 2d chain at the onset according to 



N(3) 



L 



L 



_ L p {2) 
L p (3) L p (3)L p (2) L p (3) 



L p (2) 



N(2). 



Using the dependence of the persistence length on di- 
mensionality [3g], L„(3) ~ l/(d— 1), we obtain iV(3) = 
27V(2). According to [37| N(2) m 0.75 at the onset, hence 
we obtain N(3) ~ 1.5, which is not far from our numeri- 
cal result, 2V(3) w 1.3 (see Fig. [Q. 

We now will address the issues of accuracy of the calcu- 
lations, which depends on the size of the matrices D and 
M s , and the maximal I at which the summation over I is 
stopped. In order to check our computation we verified 
that at large N (N — 100) the numerical evaluation of (|3|) 
gives with very high accuracy the Gaussian distribution 
3/(2iV) exp(— 3y 2 /4:N). The general tendency is such 



FIG. 2: Normalized distribution function G(to,x = Q,y,N) 
for chain length N = 0.5. Squares: truncation with 13 x 13 
matrices; continuous line: truncation with 20 x 20 matrices. 



that the sufficient level of matrix truncations and the 
number of terms in the expansion over the Legendre poly- 
nomials increase with decreasing N. In the limit N — > 
the whole series over I should be taken into account. The 
accuracy of the computations is demonstrated in Fig. [2] 
showing the computation of the distribution function for 
N = 0.5. The squares and the continuous curve corre- 
spond to the truncations by 13 x 13 matrices and 20 x 20 
matrices, respectively. In both cases the summation was 
stopped at I = 8. The corrections due to higher Z-s are 
negligibly small. For example, the corrections associated 
with I = 10 and I = 12 contribute only in 3rd and 5th 
decimal digits, respectively. Thus, the computations de- 
picted in Figs. [HE] can be considered as exact. 

Figure [3] shows the results of the computation of the 
3d distribution functions Go(r, N) of the free polymer ob- 
tained by performing the inverse Laplace-Fourier trans- 
form of the term I = in Eq. (TTJ) for different chain 
lengths, and its comparison with the Monte Carlo simu- 
lations [2lj| . Our results are in excellent agreement with 
the numerical data. 



IV. CONCLUSION 

To conclude, we have studied the transverse distribu- 
tion function of the free end of the three dimensional 
wormlike chain with fixed orientation and position of the 
second end using the exact solution for the Green's func- 
tion of the wormlike chain. Within the procedure of trun- 
cations of the exact formula with finite order matrices 
we find that the distribution function G(to,x — 0,y,N) 
for intermediate chain lengths, belonging to the inter- 
val 1.3L p , 3.5L p , possesses the bimodal shape with the 
maxima at a finite value of the transverse displacement, 
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which is consistent with the recent studies [13, [38| and 
[39l ] for the two dimensional chain. In contrast to the 2d 
wormlike chain, the transverse Id distribution function 
of the 3d chain shows only a tiny peak at y — in the 
vicinity of the onset of bimodality, which however disap- 
pears for larger N. We present also results of analytical 
considerations for short and large polymers which are in 
complete agreement with the classical works 0, [3, Ojj 
where these limits were investigated using different meth- 
ods. The computation of the three dimensional distribu- 
tion function of a free polymer is in excellent agreement 
with the Monte Carlo simulations plj . 



FIG. 3: The end-to-end distribution function of a free poly- 
mer. Solid lines: the distribution function of a free three 
dimensional polymer chain for N = 10, 5, 2, 1, 0.5 (from left 
to right) computed with 20 x 20 matrices. The symbols are 
the Monte Carlo simulation data extracted from Fig. 1 in 
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